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Abstract 

The BaBar experiment is a new generation detector located at the SLAC B 
factory PEP-II ring which should start taking data at the end of 1999. Its main 
goal is the study of CP violation in the B°B° system. After explaining the nature 
of this CP violation, I review the scientific program for achieving this study in 
many different modes, in the light of the recent developments obtained both on the 
experimental and theoretical side. Implications for the Standard Model are then 
discussed. 
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1 Introduction 



So far the violation of CP symmetry has just been observed in the neutral Kaon sector. 
The Standard Model can accommodate for such a violation, through the CKM mixing 
matrix. Furthermore, it even predicts CPviolation in the B°B° system. A first task of a 
B factory is thus to check whether such a prediction holds on in the B° sector. 

The CKM matrix is presently one of the less tested sector of the Standard Model. 
Indeed, two of its four parameters are presently very badly known ( p and rj in the 
Wolfenstein parametrization) . The knowledge on these two parameters is depicted in the 
so-called Unitarity Triangle (UT), where the apex of the scaled triangle is precisely the 
p,r] point (Fig |). 

n 




Figure 1: The Unitarity Triangle representation. Note that its base is normalized to 1 

Presently the three sides of the triangle are measured, but the extraction from the data 
of CKM quantities requires the knowledge of model-dependent theoretical parameters 
(coming from non-perturbative QCD and models for heavy to light transitions). Another 
constraint comes from the measurement of CP violation in the kaon system, but here 
again, due to large theoretical uncertainties, this constraint is quite weak. The use of limits 
on the B s mixing frequency is interesting but again plagued by a theoretical parameter 
(section |5T|) . 

The goal of B factories is to measure two angles of the UT (a and f3) in a clean way. 
Combining all observables will allow to (over-?)constrain the CKM matrix. Furthermore, 
time dependent CP violating asymmetries, being rare processes, are sensitive to New 
Physics phenomenon. Or, said in another way, many extensions of the Standard Model 
includes some new sources of CP violation [|IJ that could be observed at a B-factory 
experiment. 



1 



2 Which CP violation? 



In the T(4S) — > B°B° decay, after the decay of a tagging flavor B, the time distribution 
of the decay of the other B (to a final state /) is of the form: 

p{t) = Ce~ rt {l + |A/| 2 +(1 - |A/| 2 )cos(Am d t)-27mA/sm(Am d t)] (S°tag) 

p(t) = Ce~ rt [l + |A/| 2 -(1 - |A/| 2 )cos(Am d t)+27mA/sin(Am d f)] (5°tag) (1) 

where A/ = | ^[goII/) 0' ^ m d is the 7?d mixing frequency and T its width. One notices 
the difference in signs in the above expression. 

Choosing a final CP eigenstate, the time dependent asymmetry a(t) can be different 
from 0, indicating CP violation: 



a(t) = 



N(B°{t) - 


-> /) - N(B°(t) - 


-/) (i-iA/i 2 : 


)cos(Arridt) — 2Im\ fsm(Amdt) 


N(B°(t) - 




-f) 


(1 + IA/l 2 ) 



(2) 



There are two ways for the ratio to be non zero: 

• |A/| + 1 

This can be achieved either by |~| 7^ 1 or ^^H^^] 7^ 1- The former inequality 
represents a CP violation in the mixing (indirect) and the latter a CP violation 
in the decay (direct). The amount of indirect CP violation is expected to be very 
small in the B°B° system (at a level of 10 -3 ). Direct CP violation however can be 
different from in rare processes (beyond tree diagrams) and depends on the modes 
studied. 

• 7mA/ ^ 

The term 7mA/ has no reason to be equal to 0. In some "clean" cases it can 
even be directly related to the angles of the Unitarity Triangle: 7mA/ = sin 2a or 
7mA/ = sin 2(5. It arises from the interference between the decay with and without 
mixing. It is the prime motivation for the construction of B-factories. 

3 Introducing BaBar 

To achieve an experimental study of such time dependent asymmetries, the following 
requirements must be fulfilled: 

• produce a coherent 7?°7?°state (i.e. run a the T(4S)resonance). 

2 q,p appear in the physical states decomposition: \Bl) — p\B°) + q\B°) and \Bjj) — p\B°) — q\B°). 
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• since CP modes are rare (BR of the order of 10 4 , 10 5 ) have a high luminosity, 

• the time variable that appears in Eq.(^) being the decay time between the two B 
decays (t = At), it is crucial that they do not decay at the same point (otherwise the 
sin(Am^) term cannot be measured in a time dependent way and the time-integral 
over this term vanishes): one needs therefore to boost the B°B° system, i.e. use 
asymmetric beams. 

The BaBar detector is located at the PEP-II storage ring, a high luminosity collider 
(£ = 3.10 33 cm _2 s _1 is expected) of 9 GeV electrons against 3.1 GeVpositrons. This gives 
a boost to the B°B° system of •y/S = 0.56; the mean separation between the two B decays 
is about 260 /im. 

The (asymmetric) detector is a classical one for e + e~ colliders (except for the DIRC), 
made of high quality components. Going from the beam pipe (Fig 0): 

• a 5 layer silicon vertex tracker 

• a low density He-based Drift chamber 

• a CsI(Tl) calorimeter with high granularity. 

• a DIRC (Detection of Internally Reflected Cerenkov light) for particle identification 

• a superconducting coil of 1.5 T. 

• An instrumented flux return optimized for [i and K® detection. 

The DIRC is a new detector for Particle Identification based on the Cerenkov light 
emission of a particle passing through a quartz bar. While generally the light captured 
in the radiator is lost, here one uses this component which is trapped inside the quartz 
bar and propagate by internal reflection to the end, conserving its characteristic angle. 
At the end of the bar, the photon propagates into a large volume of water (the "standoff 
box") and reaches a huge array of about 13 000 photomultipliers. The reconstruction of 
the angle between the hit PMT and the bar allows to measure the Cerenkov angle, and 
thus the nature of the track. 

4 The CP program 

In order to extract CP violation parameters, one needs in real life to perform the following 
program: 
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4.1 Reconstruction of a final state / 

This is performed by usual techniques (mass peaks...), for many different final states: 

f = J/ipK° s , 7T+7T-, tt+tt^tt , 4tt, D + D~ , j/lf) K* , D* + D*~ , 1p(2S) K, .. (3) 
I will detail the first three modes while describing the following of the analysis. 

4.2 Tagging 

The goal of this part is to tag the flavor of the B meson (b or b quark?). This is generally 
performed searching for a lepton and/or kaon in the event (Fig. [3]) A sign contamination 

1- 1 + 
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Figure 3: Ways to produce a lepton and a kaon in a B° decay 

comes from secondary leptons; usually one uses a cut (as on its momentum) to enrich the 
sample in primary leptons. However in that case, one looses the information contained in 
the secondary leptons: if it is very soft, it is more likely to be a secondary lepton, so its 
sign information should be flipped. 

A tool named CORN ELIU S [IT] has been developed in the Collaboration, in order 
to combine the information of many discriminating variables associated to the lepton. 
This is achieved using various multivariate methods Q it allows to crosscheck the different 
outputs and have a grip on systematics. But much more. It allows to assign to each 
event a probability to come from a b or b quark. This probability is then input in the final 
likelihood determining the asymmetry and exploits optimally all the available information. 

The deterministic "cut" method degrades the determination on the asymmetry by a 
dilution D = e(l — 2w 2 ) where e is the tagging fraction and w the mistag fraction. Previous 

3 presently it incorporates a Likelihood analysis, a Fisher discriminant and a Neural Network 
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estimates of this quantity gave about: D = 0.33. Using the probability method allows 
to reach D = 0.41, by combining 8 discriminating variables for the lepton. 

Notice that this combination can also be used to reject the background (generally "con- 
tinuum background" ) by combining discriminating variables based on the event topology. 
CORNELIUS can provide event by event a probability to be a BB or qq event. 

4.3 Time determination 

Once a mode is reconstructed, a vertex is performed with the remaining charged tracks. 
The difference in space between both vertices represents simply / ~f/3ct where 7/3 is the 
known boost of the machine (=0.56). The resolution on the distance between both vertices 
is, for the J/ip K° s mode, about 50 jum, well beyond the mean 260 /im quoted in section 
for the mean B separation. 

4.4 Extracting CP 

There are two aspects in extracting a quantity relevant for physics. The first one is mainly 
experimental and is based on the knowledge of the detector. It will be illustrated on the 
J/ipK® mode. Going from a measured asymmetry to a relevant CKM quantity is a more 
theoretical problem, that will be illustrated on the 7r + 7r~ mode. 

4.4.1 Experimental side: introducing the /CZ7V variable (J/ipK®) 

The extraction of the asymmetry can be performed by a likelihood fit to the observed 
events. However it is more convenient to use the KTM variable. 

In a simple case (the theoretically clean mode iJ/ipK®) and neglecting for the time 
being detector effects, the event distributions (P can be written: 

B°tag:p{t) = Ce~ rt [l- sin 2/?sin(Am^)] 

B°tag:p(t) = Ce~ rt [l+ sin 2f3sm(Am d t)] (4) 

Constructing event by event the asymmetry 

Stag- — - = S t agSm(Am d t)sm2(3 = Ksin2/? (5) 
P + P 

(-) 

(where S tag = ±1 for B tag) allows to fill an histogram of the JCIAf (k) variable (see 
Fig |). 
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Figure 4: An example of the KXJ\f distribution on J/ipK® mode. From the upper box, 
one can get the whole information on the determination of sin 2(3 

The distribution of k variable has the nice property that, to a very good approximation, 
one can get the estimate of the asymmetry (sin 2(3) by the very simple formula [|J: 

This means that a single plot (as Fig ^) carries the whole asymmetry information, 
and that the asymmetry measurements can be obtained via the number of entries, the 
mean and the RMS of this histogram. Furthermore one can incorporate in the KXM , 
the tagging probability event by event and the time resolution measurement ||. All the 
JCIAf results still hold. Finally notice, that the different modes can be combined in a 
straightforward way by just summing the k histograms. 

Using the KIM approach , a recent analysis of the J/ipK®mode has been performed 
0. The measurement obtained for (one "nominal" year) is: sin 2(3 = 0.82 ± 0.15 

(while .70 was generated). 

4.4.2 Theoretical side: the penguin world ~) 

Contributing to a final state as 7r + 7r~, can exist, beside the Tree Cabbibo suppressed 
mode (a), some modes as "penguin diagrams" (b). Recently emphasis has also been put on 
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"long distance penguins" (c) (or "charming" penguins) which are of QCD non-perturbative 
nature and results from annihilation/re-scattering processes. 




Figure 5: Possible contributions to the tt+tt mode, in term of local 4-fermion operators: 
(a) tree (b)penguin (c) "charming" penguin 

The theoretical estimate of the penguin contribution is very delicate (and model- 
dependent). For this mode, it is expected that the penguin contribution is "smaller" than 
the Tree one. However, recent CLEO measurements of BR(B° — > Kir) || indicates that 
these penguin modes indeed exist and should not be neglected in extracting the CKM 
quantity sin 2a. 

Under the influence of penguins, the time distribution of the events Eq.(|I]) can be 
written as: 

Ce~ rt (l ± Rcos(Am d t) =p sin 2(a + 5)sin(Am d t)) (7) 

where R is the amount of direct CP violation, 5 is a shift due to the presence of penguins 
and a is the CKM angle. 

What one can extract from the data is a sin2a e ff but the penguin shift is unknown. 
There are several solutions to this problem, depending on what will be measured: 

• Gronau and London || have shown that measuring the decoupled amplitudes: 

A(B°^n + n-) , A(B° -> tt+tt-) 

A(B + -> TT + TT ) , A(B~ -> TT-TT°) (8) 

A(B° - ttV) , A(B°^n\°) 

allows to extract a by using the isospin symmetry. This is however performed with 
an 8-fold ambiguity on 2a. Furthermore, the amplitudes of the tt tt mode are 
expected to be small (color suppressed) and experimentally difficult to reconstruct. 

• In the case where only an upper bound has been obtained for this tt tt mode, 
Grossman and Quinn [0 have shown that the penguin shift is limited by: 

,-2 K BR(B° -> A ) + BR(B° -> tt°tt ) 

- BR(B+ -> TT+TT ) + BR(B~ -> TT-TT ) 1 ' 
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This can be very useful when constraining the penguin effects. 

• Finally one will have to rely on the theorist understanding of the penguin , reduc- 
ing the model-dependence to a minimum number of parameters |§, to obtain a 
systematic error on the determination of sin 2a. 

4.4.3 The full problem: 

A clearly challenging mode to extract a is tt + 7t~7t°. In that case the situation is compli- 
cated by the fact that: 

• p — re is assumed to dominate but p + n~ , p~n + , p°7i° interferes. 

• Experimentally the signal is not so pure (signal:background ~ 1:1) 

• There is a unknown contribution from penguins. 

However this mode is important since, it is expected to have a higher branching ratio 
than 7T + 7r~. Also since it is a non-CP final state (due to phase space) it can have a large 
cos(Am c ;t) contribution, leading to a simultaneous determination of cos 2a and sin 2a. 
This would definitely reduce the ambiguities due to a single measurement of sin 2a (in 
which case 2a = arcsin(sin 2a) and 2a = it — arcsin(sin2a) are both solutions). 

The observed asymmetry in this case is of the form: 

Ce- rt [l ± &($)cos(Am d t) =f c($)sin(Am d t)] (10) 

where &($), c($) are functions of the phase space (as the 2 Dalitz plot coordinates). 

If one collects enough data, the study of the time-dependent Dalitz plot allows to fit 
&($), c($) and extract from the data all the information on a and penguin contributions 
||. This requires however a large statistics, and as in the tt + tt~ case, the measurement 
of the color suppressed contribution (here p°7r°) is mandatory. Based on some models for 
the branching ratios || , this could require about 3 years of data taking. 

In the first year(s), the approach to this problem will be a 2 body approach: phase 
space is integrated , using relativistic Breit-Wigner, and taking into account interferences 
between p + n~ and p~n~. The effects of the penguins are neglected and will induced 
a systematic error. A recent analysis of this channel [|IU| obtains, for one year of data 
taking, an effective asymmetry: sin2a e // = 0.26 ± 0.15(stat) while the generated value 
(with penguins) was 0.43. This can give an idea of the induced penguin shift. 

Notice however that for this mode, as for 7r + 7r~ very much depends on what the 
different measured branching ratios will be. 
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4.4.4 Modes studied in BaBar 

So far, I have just described 3 analyzes. Many more channels are in fact studied and 
Table [I] summarizes the different modes which allow a determination of the angles /3 and 
a of the Unitarity Triangle. 



Angle 


Mode 


quark process 


penguins 


P 


Charmonium Kg(K®7) 
Charmonium K* 


b — > ccs 


\P\ 






D+D-, D* + D*-,DD* 


b — ► ccd 




[T] • 






b — > sss 


\P\ 


> l T l 


a 


7T + 7T~ 


b — > uud 


\p 

\T 


< 1 




7T + 7r _ 7r°(p7r) 




n? 

T\ ■ 




7r + 7r _ 7r 7r°(ai7r, pp) 






Ell 

T\- 



Table 1: Different modes studied in the Collaboration. The last column indicates the 
relative contribution between Tree processes and Penguins 



5 Implications for the Standard Model 

5.1 Present knowledge of the Unitarity Triangle 

Assuming the Cabbibo angle is known well enough, the observables that constrain the 
other 3 parameters of the CKM matrix (i.e. A,p,rj in Wolfenstein parametrization) are: 

• I Kb I which is a direct measurement of A. A lot of effort both on theoretical side 
(to understand corrections to HQET) and experimental side has been invested 
The present knowledge is (conservatively) |T2| \V c b\ = 0.039 ± 0.003 

• | y*- 1 is obtained from the exclusive or inclusive study of b — >• u decay. In these heavy 
to light transitions, the theoretical ground is much less firm. Even with a limited 
sample, the model-dependent error dominates and it is reasonable to assume for it 
a relative error as large as 25%. 

• Am<i (the mixing frequency) which is a measurement of |V^V^|. It is now well 



known thanks to the LEP time- dependent measurements [13|. In order to extract a 



relevant CKM quantity, one needs to know the theoretical parameter f Bd JB Bd |15 . 
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There exists a large spread of estimates for this value depending on the model used 
(lattices, QCD sum rules, quark models...). A reasonable range for these estimates 



is f Bi yJB Bi E [160,240]MeV 

Am s . LEP provided stringent constraints on the mixing frequency for the B s meson 
13]. To extract CKM parameters, one needs in principle to know a parameter 



analogous to the B& case: f Bs \JB Bs . However combined with Am^, knowing the 

JBs 



SU(3) flavor correction: ^ = — dV \ — = is enough. This latter is better known from 



lattices calculations. Still, recent estimates give: Q € [1-12, 1.48], and not knowing 
more than that, one must, in order to be conservative, take the upper limit of 1.48. 
Notice that LEP provided more than a limit (a set of "amplitudes" |Tj| ) and that 
this information can be used optimally, as exposed in 



• \sk\ is the measurement of indirect CP violation in the kaon sector. Here the 
QCD non-perturbative parameter Bk is quite unknown. A reasonable range is: 
B K e [0.6, 1.] 

Before combining these observables, let notice that there is a clear part of subjectivity 
for the theoretical parameters used (depending generally on personal preferences). It is 
certainly a delicate matter to estimate which model is right, and what the "error" quoted 
means? An old Bayesian ghost also appears: not knowing which model is right is not the 
same than taking a flat distribution between all estimates. 

BaBar has adopted the following way of combining , which is statistically meaningful 



The errors on a quantity are divided into 2 types: experimental errors are considered 
to be gaussianly distributed and enter a x 2 (A P> v) estimate. Theoretical parameters are 
scanned within reasonable ranges: for each scanned parameters, the x 2 is minimized 
leading to an estimate of A, p, rj and for instance a 95%CL contour in the p — r\ plane. 
A x 2 probability cut is applied in order to check the compatibility between the various 
observables. If the contour survives, one then go to the next scanned theoretical param- 
eters, etc. Knowing the exact theoretical parameter value, one could fix which contour 
is the right one. Not knowing it, one takes as a conservative choice the set of all the 
contours as the overall 95%CL knowledge of p and rj. Using this method, the present 
(1998) knowledge of p — rj is depicted on Fig. |5] (together with the values used). 

Working in another basis (A, sin 2a, sin 2(3) , the x 2 minimization can be performed 
and one obtains in the same way the overall 95% CL in the sin 2a, sin 2/5 plane (Fig. [7|). 



by "reasonable" we mean a conservative range obtained after discussion with many theorists 
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P( X 2 )> .05 




I I I I I I I I U I I I I I lUj I I I I I I I LlJ 1 I I I i I I I I I I I I I 

-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

P 



Experimental inputs 


Theoretical scan 


Am d =0.472±0.018 ps _1 


0.16<f B VB B <0.24 GeV 


l£ K l=0.002258±1.8E-05 


0.6<B K < 1 


IV ub /VJ=IV ub /VJr±0.005 


0.06<IV ub /V cb l T <0.1 


IV cb l=0.039± 0.003 


f, 2 = 1 .48 


m t =165±6 GeV 





Figure 6: The 1998 knowledge of the p — rj: one circle represents the 95% CL obtained 
when theoretical parameters are fixed. The set of all contours represents our overall 95% 
knowledge of pi] when scanning the possible range of theoretical parameters. Also drawn 
in dashed lines are the constraints provides by each individual measurement 
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1998 situation 



CO. 

'5 0.8 



0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

sin(2a) 



Figure 7: The 1998 knowledge on sin 2a — sin 2(5 using the same set of values as in Fig |] 
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From this figure, two important remarks arise: 

• sin 2/3 is presently constrained to lie somewhere on [0.4,0.8] (95% CL) []. The task of 
modes measuring sin 2(3 will therefore to test the Standard Model by checking the 
compatibility with this range. 

• sin 2a is presently unknown. The goal of a B factory is therefore to measure this 
angle. 

5.2 What a B factory can bring 

It is presently quite delicate to foresee what the impact of the B-factory will be in con- 
straining the CKM matrix, since many modes depend on what the branching ratios, 
penguins, etc.. actually are. BaBar has updated |1]J many of its analyzes with a realistic 
simulation and the full reconstruction program. A possible scenario for 30/6 -1 integrated 
luminosity (one year) giving estimated values of these angles is 

sin 2a = O.liO.l^iO.2^ (11) 
sin 2(3 = 0.6±0.08 ea; p (12) 

These numbers represent a reasonable order of magnitude but the central values are 
completely hypothetical (within the Standard Model). Fig[| shows the impact of such a 
measurement in the p — t] plane. 

Two final remarks: 

• For this combination one neglects the improvement that will appear with time on the 
present observables. In particular, \V c b\, |y |, Am^ measurements should improve 
with B-factories. 

• Just as knowing the sides of the UT allows to already constrain one of the angle (Fig. 
|7|), in the other way, measuring sin 2a and sin 2(3 will allow to measure in a model in- 
dependent way, the theoretical parameters (fB d \jB Bd , B K ,\y^\ t h). This happens on 
Fig. H where the combined contour is not better than the simple BaBar constraint. 

6 Conclusions 

Several tools have been developed in the BaBar collaboration these last years. In partic- 
ular: 

5 Recall however than this is not a measurement and relies on a set of preferred theoretical parameters: 
no density distribution can be derived for it (since the distribution of theoretical parameters is unknown) , 
just a range. 
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Experimental inputs 


Theoretical scan 


Am d =0.472±0.018 ps" 1 


0.16<f s VB 8 <0.24 GeV 


l£ K l=0.002258±1.8E-05 


0.6<B K < 1 


IV ub /VJ=IV ub /V cb l T ±0.005 


0.06<IV ub /V cb l T <0.1 


IV cb l=0.039± 0.003 


^ = 1 -48 


m t =165±6 GeV 




sin 2a= sin(2a) T ±0.1 


-0.1<(sin2a) T <0.3 


sin 20=O.6±O.O8 





Figure 8: A possible scenario for the determination of p — rj after one year of BaBar data 
taking. Also shown as a new contour, in light grey (or blue, in color) is the BABAR-alone 
95% CL. 
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• The tagging is now a probabilistic answer of a multivariate analysis based on dis- 
criminating variables. This incorporates optimally the full information extracted 
from an event. 

• The KXM variable is a golden one for the extraction of an asymmetry from the 
data, modelizing all detector effects. It allows to present results in a clear way, to 
optimize selections and combine the results of many channels (Collaborations?). 

• A method has been developed to combine observables constraining CKM parameters 
in a statistical meaningful way. 

A B-factory will provide precision tests of the Standard Model. First of all, it should 
prove that CP violation exists in the i?°i? sector. Then many channels will be combined; 
this require a deep understanding of the detector. Then, in conjunction with theoretical 
work, the relevant CKM angles will be extracted. 

Presently sin 2/3 is bound to lie in [0.4, 0.8] (@ 95%CL) from indirect measurements 
(this however implies the choice of theoretical model-dependent parameters). Therefore 
the goal of a B factory for this angle is to check whether the time- dependent asymmetry 
is compatible with this range. The golden modes are Charmonium K° s and Charmonium 
K* . Since these modes are theoretically well under control, any significant deviation from 
this range would indicate a problem (on errors? theoretical parameters? New Physics?). 
In particular, if no asymmetry is observed in J/ipK° s mode, this would rule out the CKM 
mechanism of mixing between 3 fermion families and indicate New Physics. 

Presently sin 2a is unknown and its determination is a challenging task for B-factories. 
Very much depends on the values of the branching ratios. In particular if B° — > 7r°7r° 
(resp. B° — > p°7r°) is measured it allows a model-independent determination of sin 2a 
from 7r + 7r~ (resp. 7r + 7r~7r° ). The many accessible rare modes which will be studied at 
B-factories will allow to test different models (as factorization or SU(3) symmetry from 
B(B° — > Kn)...) and get an insight into the unknown world of penguins. 
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